UNIVERSITY OF ILLINOIS 
DIGITAL COMPUTER 



TITLE . 

TYPE 

METHOD OF USE 



CAPACITY 



DURATION 

PUNCHING OF THE TAPES 



LIBRARY ROUTINE K Ik - 233 



By Gene H. Golub 
Multiple Regression Analysis with Transformations 

Entire Program 

1) Read in master program. A sum check is performed, 
and if the program is read correctly the machine will 
stop on 2k (012)^; otherwise it will stop on an FF order. 

2) Read in the parameter tape. The computer stops on 
20 (0NL) l6 . 

3) Read in the data tape. 

k) After all data have been read into the computer, the 

desired information will be printed out. The computer 

stops on 2k (012 ) /- . 

5) Another problem^ can be begun by reading in new 

parameters-. 

The maximum number of nonr: eliminated variables is twenty-two. 

The total number of variables must be < 50. The number of 

observation is unlimited.. 

This program hangs up unless s > n + 1 * 

[See Note 3] 

For each problem six parameters are necessary. 

They are as follows: 

(1) Let "s" be the sample size. Put sS on the tape. 

(2) Let "n" be the number of non-eliminated variables. 
Put nN on the tape. 

(3) Let "f " be the number of decimal places to which 
correlations are to be printed . Put fF on the tape . 
(k) Let " JL" be the number of decimal places to whicn 
covariances are to'be printed. Put J2L on the tape. 
(5) Let "k" be the number of decimal places to which 
standardized regression weights, their standard errors, 
and multiple correlations are to be printed . Put kK 
on the tape. 
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(6) Let "J* be the number of decimal places to which 
unstandardized regressions weights, their standard errors 
and standard errors of estimate are to be printed. Put jJ 
on the tape. 

Any of the four parameters f , j, k, or £ may be set ±6 aero. 

erea ses- 
e 1 the first variable will 



The parameters are followed %'a sequence of'l'a ami 
l*s indicate vhich of the non-eliminated variables are to be cdnsTd 



glvely independent, (e.g., if this sequence is a si _ 

be considered dependent, if it consist of 1001 the first variable <4nd then the fourth 
variable will be considered dependent), d's following . the last 1 in the sequence 
should be omitted. If more than one variable is to be considered dependent no delays, 
spaces, line feeds or other characters withHhe fifth hole punched may be inserted 
between the first 1 and the last 1 in this sequence.? 8 The first character after the 
last 1 in tile sequence must be a character with the fifth' hole punched^ (e.g., delay, 
space, line feed, etc). The portion of the input tape between this character and 
a character which consists of a single fifth hole on the tape will be exactly reprcj* 
duced on the output tape. In this interval some identification may be punched. The 
identification can be a number, a name, or any combination of characters which does 
not include a character consisting a single fifth hole on the tape. It must, however, 
be followed by such a character* If no identification is desired, two fifth hole 
characters most be punched immediately after the last 1 in the sequence. The second 
of these must consist of a single fifth hole on the tape* Ho fifth hole character 
say be punched following the S, V, F,'L, J, K parameters in the section where the 
successive dependent variables are Indicated. If one is punched, characters between 
the fifth hole character the single fifth hole delay will be reproduced as identifica- 
tion. 

Next, a sequence of signed integers, followed by an H, is 
punched on the tape. These integers take on the values 0, 1, 2, 3, k 9 5, 6. 
Corresponding to each integer, the following transformations are made: 



n 


transformation on x 


+0 


eliminate variable 


+1 


X 


+2 


x 2 


+3 


X 5 


44 


'VxT x > 


+5 


log x/10, x > 10" 10 


46 


-- arc sin x 

JC 
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The number of digits in this sequence must be equal to 
the total number of variables. 

Each element of an observation vector must be punched 
as a sign followed by up to 12 decimal digits. The character N must be punched 
af+er each row. If an F follows a row, then the machine will stop, and another 
part of the data tape may be placed in the reader. 

SCALING OF THE DATA 

Each observation must be scaled so that it is less than 
one in absolute value before it can be punched. Each element should be scaled 
"by the same value if the constant term is to be easily interpreted. If each ele- 
ment is scaled by a constant c, then the correlation coefficients, multiple 
correlation, regression coefficients and standard error of the regression co- 
efficients will be unaffected.' The means, standard deviations, standard estimate 

of error and constant term will all he scaled by c. The variances and covariances 

2 
will he scaled by c * 

THE PRINT OUT 

This is best illustrated by a typical example For 
a problem with four variables, the print out could 
appear as follows on the next |?age s The means and 
standard deviations are always printed to nine places . 
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11- 



'21- 



'22' 



Correlations 



Means 




+a.ooooooo +.01386883 



+0.l64934l 
+1.0000000 

+0.83^721 
+0.5308620 

+0.3013304 

-0.4067406 

+0.2975063 
+1.0000000 

+.214150000 
+.549600000 
+.106490000 
+.105000000 



+.00459466 
+.05595594 

+.00510478 
+.00652617 

+ ,00270089 

+.00204625 
-.00554800 

+.00089155 
+.00332500 

+.117765986 
+.236550079 
+.051970115 
+.057662813 



Covariances 



■Standard deviations 



Dependent 
Variable 



Constant term. 



+.043813563 
+.928217019 

-1.00000000 
-0.33097420 
+2.92808434 
-O.72196344 
■^ |+0.l6004788 1 

Unstandardized 



■Standard error of estimate.. 
■Multiple correlation 



+.000000000 
+.076293334 
+.332283964 

r + .2778153451 
? 

Standard 



-1.00000000 

-0.66480972 

+I.29216327 

1-0.353501411 



+.000000000 

+.153246236 
+.146636873 

I +.156029212 



Regression Weights Errors 



Standardized ' Standard 
Regressions Weights Errors 



SIM CHECK 



Throughout this program, there are a number of checks 
on the tapes and magnetic drum storage. They are as 
follows : 

1.) If the master tape has been read in incorrectly, 

the machine will stop on an FF Oil order at (378)^ 

2 • ) If the data tape does not have the same number 

of variables in each row, the machine stops on an 
FF 012 order at (0j4) ^. 



MATHEMATICAL MODEL 
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3») If a drum error has been made, the machine will 
stop on an FF 013 order at (o6x) s. 

Let x be the jth observation on the ith variable. 
As each observation vector is read in, 

I jzhL and l !m k<i 

J=1 s j=1 s ill, a, ...n 

k = 1, . . ., -n '' ■ ** 

are computed. After all the observations have been 

read in, the correlation matrix is calculated and then 

inverted using a modified version of M Ik. The formulas 

used for calculation are as follows: 

s x,. . 
The mean: x = Z — ii 

j=l s 
The Covariance: s x x 

ik A s - x * ** 

The Correlation Coefficient: c 

ik 



r ik = 



Jc, c 



ii kk 



The Standardized Regression Weight: 

. ik 
^ = i| * k = 1, ..., n 
r 

ik 

where r is the (i,k) element of the inverse of the 

correlation matrix. 

The Standard Error of the Standardized Regression Weights 



\ 



* 
ik 



NJ 



r ii r a , (r ik )2 

(s - n) (r 11 ) 2 
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The Multiple Correlation Coefficient 



R. = 

1 , 



1 - 



1 



li 



NOTE 1 



NOTE 2 



NOTE 3 



Number 
of Variable; 



10 
15 



The Standard Error of Estimate! 



S. =|c. (1 - R, ) 
l J n l 



The Unstandardized. Regression Weights 



lk 



N 



li 



'kk 



b* 

lk 



The Standard Error of the Unstandardized Regression 
Weights: 

S, 



b 



ik N 



ii 



'kk 



ik 



Beginning at drum location 25o0, the subroutines 
necessary for this program are stored. 
Beginning at drum location 358 ] +, the means, standard 
deviations, correlation matrix and its inverse are 
stored . 



A TABLE OF TYPICAL TIMES 



Read in 100 
Observations 
of 8 digits 



45 sec. 

95 sec. 
150 sec . 
218 sec. 



Print Correlations 
and Covariance 
Matrix to 8 Places 



10 sec 

32 sec 

68 sec 

120 sec 



In addition, it requires about 35 



Invert 

Correlation 

Matrix 



Print Standardized 
and Unstandardized 
Regression Weights 
and their Standard 
Error to 8 Places 



k sec . 
10 sec. 
18 sec . 
35 sec . 



6 sec o 
12 sec . 
18 sec. 
23 sec. 



■econds to read th e marter tape into the machine 

Kt.s 7/1 /fin H 



DATE March 20, 1957 
PROGRAMMED B Y y<5 /&rt~J^ 
APPROVED m /jZ, 7yiu£fjtA_ 



LOCATION 


ORDER 


NOTES PAGE 1 K Ik 


Decimal Order Inpu 


t- X-l 






00 3K 






3 


00 F 




Location of data, means and matrix 




00 500F 




of cross-products 


k 


00 F 

00 960F 




Location of Library Routine N 12 


5 


00 F 

00 58F 




Location of List of Constants 


6 


00 F 
00 ll6F 




Location of Library Routine R-l 


7 


00 F 








■: 00 358^F 




First location of Data Stored on Drum 


8 


00 F 




Location of Diagonal Elements of Inverse 




00 982F 




and Scaling Factors . 


9 


00 F 

00 76F 




Location of Library Routine Y-l 


10 


00 F 
00 12£F 




Location of Library Routine P-l (modified) 


11 


00 F 








00 2560F 




First Location of Subroutines on Drum 


12 


00 F 








00 700F 






13 


00 F 








00 IfOOF 




Location of directory 


14 


00 F 








00 275F 




Location of transformation routine 




00 76K 






Library 
Subrout 


Routine Y-l 
ine A 




Transfer Blocks of Words from the Memory 
to the Drum or from the Drum to the Memory 
Change FF order in Y-l here to read FF 19F 


This subroutine clears the memory, inputs the parameters, and translates the 


identification . 








00 200K 









K5 F 






. 


k2 40L 




Plant link 



LOCATION 


ORDER 




NOTES PAGE 2 




1 


L5 *OL 
to 2L 








2 


00 11F 


from h 








kl ( )F 


by 1,3 


Clear memory 




3 


F5 2L 
to 2L 








k 


LO lUL 
36 2L 








5 


92 135F 
41 7F 


from 13 






6 


81 ^F 
LO 13S5 


from 10 






7 


32 10L 
Lk 13S5 




Read in parameters and 




8 


50 7F 
7^ 13S5 




convert 


J 


9 


S5 F 








* 


to 7F 








10 


26 6l 
k2 11L 


from 7 






11 


L5 7F 
to ( )f 


by 10 






12 


F5 8f 
to 8F 








13 


LO ^2L 
32 5L 








Ik 


to 8f 

81 4f 








15 


FO 5S5 
32 17L 


from 19 


Read in sequence of l's and O's 




16 


L5 9F 
l£ 8f 




to determine dependent variables 




17 


to 9F 
L5 8F 


from 15 






18 


10 IF 

to 8f 




— » ■ — ■ ' 





K Ik 



LOCATION 



ORDER 



NOTES 



PAGE 3 



K l4 



19 
20 
21 
22 

2k 

25 
26 

27 
28 

29 
30 
31 
32 
33 
3k 
35 



91 kF 
36 15L 
26 25L 

L5 7F 
32 22L 
50 24L 
26 23L 
50 23L 
04 6f 
42 24L 
02 IF 

92 ( )F 

91 kF 
kO 7F 

L7 7F 
32 20L 

92 707F 
92 131F 
92 515F 
L5 2F 
42 0S5 
k2 6S5 
42 16S5 
00 20F 
46 S5 
L5 S5 
L4 1S5 
40 2S5 
Lk S5 
40 3S5 
L5 4F 
42 8S5 
L5 5F 
42 9S5 



from 26 



from 21 



by 23 



from 20 



read in identification and print out 



Set parameters 



LOCATION 
36 

37 
38 

39 
kO 
kl 
k2 
kl 



r 



ORDER 

L5 OF 
1*2 10S5 
L5 3F 
1*2 11S5 
L5 IF 
kO 15S5 
L5 9F 
kO 1^S5 
00 IF 
22 ( )F 
80 11F 
1*1 103^F 
80 F 
00 6f 
00 11F 
kl S3 



NOTES 



PAGE k 



K lk 



by 



Store subroutine A 
on drum 

[This is an interlude which is 
later overwritten.] 



00 800K 
J JO 200F 
50 L 
26 S9 
00 SS 
00 kkY 
L5 OF 
kO kL 
26 999F 

26 800N 
Subroutine B 
This subroutine computes the means and sum of cross-products divided by s 



00 200K 
S5 50F 
kO OL 
52 SJ 
50 1L 
26 Sk 
L5 21S** 



Read in integers determining 
transformations 



LOCATION 


ORDER 


NOTES PAGE 5 K 14 


3 


4o 56L 

LO 6L 






4 


10 20F 
42 53SF 






5 


50 L 
24 7L 






6 


40 SJ 
L5 21S4 






7 


50 8L 
26 999F 






8 


00 F 
00 7L 






9 


00 F 
26 TL 
26 IN 









L5 15S5 
L4 5S5 






1 


40 4F 








K5 F 




Plant link 


2 


42 48L 
19 38F 






3 


66 15S5 








S5 F 




l/s 


4 


40 5F 
L5 49L 






5 


4o 6f 
l4 7L 






6 


46 6f 
41 9F 






7 


50 49F 


from 39 


Read in and transform 




50 TL 




observation row 


8 


2,6 SF 
LO L2S5 






9 


30 9L 








L5 21S4 




1 



LOCATION 



10 
11 
12 

13 
14 

15 
16 

17 
18 

19 

21 
22 

■ ; -23 
24 

26 



ORDER 



LO 6F 
40 OF 
L3 OF 
32 12L 
FF 18F 

L5 1S5 
42 15L 
L5 2S5 
42 16L 
46 16L . 
50 5F 
7J ( )F 
L4 ( )F 

^P ( )F 
F5 15L 
42 15L 
L5 16L 
L4 4S5 
40 16L 

F5 OF 
40 OF 

LO 6S5 

36 15L 

F5 5S5 

40 2F 
L5 1S5 

42 ZGL 

<*, 

L5 3S5 
42 29L 
46 29L 

L5 1S5 
42 28L 

41 IF 
50 ( )F 



NOTES 



PAGE 6 



K 14 



from 21 

by 13,17 
by 14,19 

by 1^,19 



Have n observations been 
read in? 



x 



ik 



k=l 



P = 1, 2, 



from 36 



by 23,35;from 3^ 



LOCATION 



ORDER 



NOTES 



PAGE 7 



K l4 



27 
28 

29 
30 
31 
32 
33 
34 
35 
36 

37 
38 
39 



ill 
k2 
43 
kk 



7J 5F 

40 OF 
50 OF 

7J ( )F 
L4 ( )F 

ko ( )f 

F5 28L 
k2 28L 
L5 29L 
L4 4S5 
40 29L 
F5 IF 
kO IF 
LO 2F 
32 26L 
F5 26L 
42 26L 
F5 2F 
40 2F 
L5 IF 
LO 6S5 
36 25L 
F5 9F 

40 9F 
LO 4F 
36 7L 

41 OF 

L5 3S5 
k2 42L 

L5 1S5 
k2 k^L 
L5 ( )F 
22 k3L 
kO ( )F 
F5 43L 
k2 kjL 



by 25,30 
by 24,32 
by 24,32 



X ik x Jk 



k=l 



p = 1, 2, 



Have s observations been read in? 



by 4l,46;from 48 
by 42,44 



Store diagonal elements 
consecutively 



LOCATION 


ORDER 




NOTES PAGE b K 14 


4 5 


F5 ^+2L 
F4 OF 






46 


42 42L 
F5 OF 






hi 


40 OF 
LO 6S5 






48 


32 42L 








22 ( )F 


by 2 




4 9 


40 S3 
L5 21S4 






50 


00 F 
00 52L 






51 


50 50L 
26 999F 






52 


00 IF 
26 51L 
26 IN 






Subroutine C 






This subroutine compute 


js the standard deviations 





K5 F 








42 l4L 




Plant link 


1 


L5 2S5 
42 5L 






2 


46 5L 
L5 1S5 






5 


42 8L 
46 6l 






4 


4l OF 

4l 3F 






5 


50 ( )F 


by 1,10? from 13 





79 ( )F 

L4 ( )F 


by 1,10 
by 3,H 


s of \o 

2j .. x.. 1 Zj x., I 

k=l ik \ ik / 

_. s 'ik _ \^s / 




32 TL 




V 



LOCATION 


ORDER 


NOTES PAGE 9 K 14 


7 


23 8L 




If the variance is negative, replace 




50 7L 


from 6 


by zero* 


8 


26 s6 




Compute standard deviation 


9 


40 ( )F 

L5 5L 
L5 4S5 


by 3,11; 
from 7 




10 


40 5L 
LO S5 






11 


46 6l 
42 8L 






12 


F5 3F 
40 3F 






13 


LO 6S5 
36 5L 






14 


00 IF 








22 ( )F 


by 




Routine H 








Read in data and perform transformations 




00 275K 









50 50SJ 
50 OL 






1 


26 S4 
40 13F 






2 


4l 10F 
L5 17L 




Read in data 


3 


42 9L 
L5 49L 






4 


42 7L 
00 IF 






5 


L5 21L 
42 6L 






6 


L5 16L 


from l4 






l4 sj 


"by -5,12 




7 


42 8L 








L5 50SJ 


by 4,11 _ 


1 , , _ _ , , . , _^T 



LOCATION 
8 

9 
10 
11 
12 

13 
Ik 

15 

16 

17 
18 

19 

20 
21 
22 
23 
2k 

25 



ORDER 

kO IP 

26 ( )F 

7J IF 

**0 S3 

F5 9L 

k2 9L 

F5 7L 

**2 7L 

F5 6L 

te 6l 
F5 1QF 
kO 10F 
LO 53L 
36 6L . 
L5 13F- 
22 215F 
26 11L 
00 16L 
22 9L 
00 S3 
50 IF 
26 9L 
50 IF 
26 kn 

kl OF 
22 k?L 
22 kjL 
00 SJ 
kO 3F 
50 3F 
7J 3F 
kO 12F 
iA- 5kL 
kO IF 
kl OF 
U IF 



by 7 

from Q,k2 
by 3,10;from rf, 1 *, 1 ^, 35,36 



from 16 



HOTBS 



Determine transformation 



PAGE 10 



K Ik 



Jump to subroutine B 
eliminate x 



log 10 x/10 



— arc sin x 
it 



x 2 ->12 
80 FOOF 



Locmoff 


ORDER 


NOTES 


PAGE 11 


26 


.32 26L 


,•' 






27 


50 26L 
26 S6 






1 


Jl-x 2 




50 5^L 








28 


LJ 12F 
36 37L 






1 

! 


29 


L5 2F 
66 3F 






i 


50 


S7 F 
50 30L 




arc tan 




31 


26 366f< 
kO OF 








32 


50 OF 
79 51L 




V* 




35 


lb 501. 
kO IF 




1/2 




3* 


L5 3F 
36 36L 


* 






35 


LI IF 
22 % 








36 


L5 IF 
22 9L 


from 3k 






37 


L5 3F 
66 27 


from 28 






38 


35 F 

50 38L 




> 




39 


26 566F 
40 IF 




arc tan 




ko 


50 51L 




I 




kl 


26 9L 
7J IF 
kO OF 


from 19 


> 




k2 


50 OF 
26 9h 




- x' 


, 


kl 


00 IF 






j 




50 k?L 


from 20 







K Ik 



LOCATION 

*T 
h8 

±9 

50 

51 

52 

53 

5* 



ORDER 

22 S6 
22 9L 
52 ^5L 

50 45L 

26 330F 

40 OF 

19 5F 

50 OF 

7^ 52L 

00 5F 

22 9L 

00 50SJ 

40 F 

00 F 

00 F 

00 3183 0988 6l8^J 

06 F 

00 okjk 29^ 8195J 

80 F 
00 ( )F 
80 F 
00 F 



NOTES 



PAGE 12 



K Ik 



lp g e X 
52 



lo log io 



1 



"by word h 9 Routine B 



00 330K 
Routine S-5 - l/32 Natural Logarithm 

00 366K 

Routine T-k - New arc Tan 

00 810K 
JO 200F 
50 L 
26 S9 
00 700SS 
00 191F 
L5 OF 
kO kL 
26 999F 



Store subroutines 
B, C, and H on drum. 

[This is an interlude which is 
later overwritten.] 



LOCATION 



ORDER 



26 810N 



NOTES 



PAGE 13 



K ik 



Subroutine D 

This subroutine computes and prints' 



correlations, covariances means, and 



standard 


deviations . 
00 200K 









K5 29L 








42 56L 




Plant lin 


1 


L3 8S5 








36 6l 




Should co 


2 


50 8S5 
75 13S5 






3 


K5 F 

iA 4s5 






k 


00 20F 
k6 28l 






5 


L5 26l 
22 6L 






6 


L5 OL 


from 1 






46 27L 


from 5 




7 


L3 9S5 
32 11L 




Should cc 


8 


50 9S5 
75 13S5 






9 


K5 F 

00 20F 






10 


k6 31L 
L5 32L 






11 


26 12L 








L5 59L 


from 7 




12 


k2 29L 
F5 5S5 


from 11 




13 


kO 5F 
L5 3S5 






ik 


42 20L 
k2 26L 







LOCATION 


ORDER 




NOTES PAGE 14 K 14 : 


15 


L5 2S5 
46 20L 






16 


L5 1S5 
46 22L 






17 


4l 6f 
L5 2S5 


from 44 




18 


42 19L 
L5 1S5 






19 


42 21L 








X5 (0)F 


by 18^34| from 39 


20 


.79 t )F 


by 15, 42 






L4 ( )F 


by 14, 37 




21 


40 3F 








50 ( )F 


by 19, 35 




22 


7.1 ( YE 1 

40 OF 


by 16, 41 




23 


L3 OF 
32 27L 






24 


50 3F 








75 57L 




Compute correlations coefficient 


25 


66 OF 
S5 F 






26 


50 28L 








40 ( )F 


by 14, 36 




27 


22 ( )F 


by 6 






23 26L 


from 23 




28 


52 ( )F 


by 4 






50 28L 




Print correlations 


29 


26 SK 








26 ( )F 


by 12 




30 


92 963F 
L5 3F 






31 


52 ( )F 
50 31L 


by 10 


Print eovariancesn 


32 


26 SK 








50 5QL , 




Waste 



LOCATION i 


ORDER 


NOTES PAGE 15 K 14 


35 


92 151F 
92 515F 






54 


F5 19L 
42 19L 






55 


F5 21L 
42 21L 






56 


F5 26L 
42 26L 






37 


42 20L 

F5 6f 






38 


4o 6f 

LO 5F 






39 


52 19L 
92 151F 






40 


L5 22L 
L4 4S5 






41 


46 22L 
L4 S5 






42 


46 20L 
F5 5F 






45 


4o 5F 

L5 6f 






44 


LO 6S5 
36 17L 






*5 


L5 2S5 

42 4tl 






46 


L5 1S5. 

42 49L 






47 


4l 5F 








L5 ( )F 


"by ^5 5 52 1 from 56 


48 


52 91F 
50 48L 




Print means 


49 


26 SK 








L5 ( )F 


by 46 5 55 . 


. 



LOCATION 


\ 

ORDER 


NOTES PAGE l6 


50 


52 91F 

50 50L 






51 


26 SK 
F5 ^7L 




Print standard deviations 


52 


k-2 kfL 
F5 ^9L 






55 


k2 ^9L 
92 131F 






5^ 92 515F 








F5 3F 






55 


40 3F 
LO 6S5 






56 


32 kjL 








22 ( )F 


by 




57 


00 F 








00 1000 000 0000J 




58 


L3 8S5 
36 3^L 






59 


26 33L 

' 00 58L 








00 820K 




■ 





JO 200F 
50 L 






1 


26 S9 
00 1?9SS 




Store subroutine 

- 


2 


00 60F 
L5 OF 




b on drum* 


5 


ko kL 




[This is an interlude which is 




26 999F 




later overwritten . ] 




26 820N 






i 

Subroutine E 






This subroutine stores the standarc 


1 deviations, means and entire 


correlation matrix on the drum. 

1 1 





K 1 



LOCATION 


ORDER 


NOTES 


PAGE 17 




00 200K 











K5 F 
42 28L 




Plant link 




1 


L5 S5 
46 22L 








2 


L4 S5 
46 7L 








3 


L4 4S5 
10 20F 








4 


l4 6l 

40 21L 








5 


JO S3 
50 5L 




Store standard deviations 
means on drum 


and 


6 


26 S9 
00 S7 








7 


00 ( )F 
L5 3S5 


by 2 






*8 


4o 6f 

42 11L 








9 


F5 5S5 
4o 5F 








10 


L5 1S5 
42 12L 


from 






11 


41 3F 










L5 ( )F 


by 8,l4,17,24;from 19 




12 


32 12L 
40 ( )F 


by 10, 13 






13 


F5 12L 

42 12L 




Square matrix 




14 


F5 11L 
42 11L 








15 


F5 3F 
LO 5F 








? 16 

• 


36 18L 
L5 3F 









K 14 



LOCATION 



17 
18 

19 
20 
21 
22 

23 
24 

25 
26 

27 
28 



ORDER 



L4 11L 
42 11L 
F5 3F 
40 3F 
LO 6S5 
32 OIL 
JO S3 
50 20L 
26 S9 
00 ( )F 
00 ( )F 
L5 6F 
L4 5F 

.4o 6f 

42 11L 
F5 21L 
L4 3F 
40 21L 
F5 5F 
40 5F 
l4 5S5 
FO 6S5 
36 10L 
22 ( )F 

00 830K 

JO 200F 
50 L 
26 S9 
00 240SS 
00 29F 
L5 OF 
40 4l 
26 999F 

26 83ON 



NOTES 



PAGE 18 



K 14 



from 16 



by 4, 25 
by 1 



Store successive rows of correlation 
matrix on drum 



by 



Store subroutine E on drum 

[This is an interlude which is 
later overwritten.] 



LOCATION 



ORDER 



NOTES 



PAGE 19 



K lk 



Subroutine F 

This" subroutine reads in rows of tlie correlation majtrix off the drum, 

computes the inyjers'e and stores the rows'of the 'inverse matrix on theudrum. 

00 ii6k 







8 
9 

10 

11 

12 

13 
14 



K5 F 


k2 l6L 


ki 8f 


L5 


S5 


k2 


15L 


k6 i4l 


k6 


1% 


hk 4S5 


k6 26l 


Lk 


S5 


k2 8f 


L5 


32L 


IX 8f 


ko 18L 


50 16S5 


L5 


1*5 


7k '16S5 


S^tf 


L4 i8l 


ko 


25L 


ko 8f 


L5 


>3k 


ko 


2g& 


F5 


33L 


14.2 


23L 


LV 


S^ ? 


46 


23L 


00 


if ; 


JO 


( )F 


L5 


jM 



by 2 



Plant link 



LOCATION 
15 

16 

17 
18 

19 
20 
21 
22 
23 

25 
26 

27 
28 
29 
30 
31 



ORDER 

26 35L 
00 ( )F 

00 IF 
22 ( )F 
50 150L 
50 17L 
26 S9 
00 F 
00 ( )F 
F5 18L 
Lk 16S5 
1*0 18L 
22 56L 
00 IF 
L5 ( )F 
1*0 ( )F 
L5 ( )F 
1*0 ( )F 
JO 150L 
50 2l*L 
26 S9 
00 F 
00 ( )F 
F5 23L 
1*2 22L 
F5 22L 
1*2 23L 
Ll* 1*S5 
k6 22L 
F5 2% 
Fl* 16S5 
1*0 25L 
26 ll*lL 
00 F 



NOTES 



PAGE 20 



K ll* 



by 2 
by 



by 6,20 
by 6,20 
by 3 



by ,^1,29 
by 11,27 
by 13 
by 12,28 



by 9,30 
by 9,30 
by 1* 



Read in successive rows of 
correlation matrix from drum. 



Store successive rows of inverted 
corf elation matrix on drum. 



LOCATION 


ORDER 


NOTES PAGE 21 K 14 


32 


26 S9 
00 S? 






33 


L5 l^OL 

4o s8 






34 


00 F 
00 35L 






35 


50 34L 
26 999F 






36 


00 F 

26 35L 
26 IN 






Routine 


'm-14 (modified) 




[Note: M-i4 has been modified so that 




00 840K 




the columns of the inverted matrix are 





JO iloF 




scaled by powers of 2 rather than by 




50 L 




powers of 10 0] 


1 


26 S9 




Store subroutine F on drum 




00 270SS 




[This is an interlude which is later 


2 


00 150F 
L5 OF 




overwritten 0] 


3 


4o 4l 
26 999F 
26 840N 






Subrout 


ine Go 




„ 4- -f-"W~ ^A-^-^A^-^A ~-~-».~-», ^v-P -»^r™«^o-! ei-r* +Viq 



This program computes and prints out the standard error of regression, the 
multiple correlation^ the standardized and unstandardized regression -weights 
and their standard errors 



00 200K 
K5 ( )F 
42 154L 
50 10S5 
75 13S5 



by 3 



Plant link 



: r 

LOCATION 


ORDER 


NOTES PAGE 22 


2 


K5 ( )F 




Compute print parameters for 




00 20F 




standardized regression weights and 


3 


46 OL 
46 123L 




standard errors. 


4 


L4 4S5 
46 2L 






5 


50 11S5 




Compute print parameters for 




75 13S5 




unstandardized regression weights 


6 


K5 ( )F 

00 20F 


by 7 


and standard errors . 


7 


46 6l 
46 119L 






8 


l4 4s5 
46 12L 






9 


00 IF 
L5 15S5 






10 


LO 16S5 
40 OF 






11 


19 38F 
66 OF 






12 ; 


S5 ( )F 

4o i4f 


by 8 


l/s-n 


13 


L5 S5 
L4 S5 






l4 


46 19L 
LO S5 






15 


l4 4s5 
46 31L 






16 


L5 8f 
40 30L 






IT 


50 SN 




Read means and standard deviations 




50 17L 




from drum. 


18 


26 S9 
00 S7 







K 14 



j LOCATICK 


ORDER 




NOTES PAGE 23 


K 14 


19 


00 ( )F 
L5 17S5 


by 14 






20 


42 23L 
F5 17S5 








21 


4^25L 
L5 56L 








22 


42 24L 










L5 14S5 


from 155 


Dependent variable? 




23 


36 148L 










L5 (t i± )F 


by 20,148 


Diagonal element of inverse 




24 


4o 9F 










L5 fe i )F 


by 22,150 


Standard deviation 




25 


40 IGF 










V? Ce ± )P 


by 21,149 


Scaling factor 




26 


40 11F 
L5 156L 




1/10 




27 


4o 7F 
L5 5S5 




-1 




28 


40 8F 
32 29L 




Waste 




29 


50 50SN 
50 29L 




Read first row of inverse matrix 
from drum 




30 


26 S9 
00 S7 


by 16,152 






31 


00 ( )F 
L5 OL 


by 15 


Set print parameters 




32 


46 I36L 
L5 2L 








33 


46 134L 
L5 6L 








3^ 


46 13OL 
L5 12L 








35 


46 128L 
46 i46l 








36 


4l 15F 
L5 56L 


from 162,166 







LOCATION 
37 

38 

39 
to 
4l 
42 

43 
44 

45 
46 

4 7 
48 

49 
50 
51 
52 
53 
5^ 



ORDEK : 

42 49L 

L4 S5 
42 85L 
L5 168L 
42 48L 
l4 155L 
42 84l 
42 127L 
L4 155L 
42 95L 
42 129L 
L4 155L 
42 57L 
42 133L 
L4 155L 
42 77L 
42 13% 

L5 1TS5 
42 50L 
F5 17S5 
42 51L 
4l 17F 
4i i6f 
L5 (t. .)F 
40 3F 
L5 (s )F 
4o 4f 
L5 (t )F 
4o 5F 

L5 (c )F 

4o 6f 

50 7F 

71 3F 

40 OF 

L3 9F 

L6 OF 



NOTES 



PAGE 24 



K 14 



Waste 



by 39,96;from 106 
*y 37,97 
ty 46,98 
fcy 47,99 



LOCATION 


ORDER 


NOTES PAGE 25 1 


55 


36 158L 




Change scale 




L5 OF 




Compute standardized regression weight 


56 


66 9F 
S5 SN 




and store 


57 


Jj-0 12F 
I4-O (b* )F 

!<5 






58 


L5 ^F 
LO lOF 






59 


kO OF 
L3 OF 






60 


32 77L 
50 3F' 






61 


79 3F 
40 OF 






62 


50 OF 








7J 6f 




Compute standard error of standardized 


63 


l*-0 OF 
50 5F 




regression weights 


6k 


7J 9F 
kO IF 






65 


50 IF 
7J 11F 






66 


lA OF 
40 OF 






67 


50 OF 
7J l^F 






68 


40 OF 
50 9F 






69 


7J 9F 
kO IF 






70 


50 IF 
7J 6f 






71 


kO IF 
50 OF 






72 


71 8F 

kO 2F 







K ik 



LOCATION 


ORDER 


NOTES page 26 


- 


73 


L7 2F 
L2 IF 








7^ 


32 163L 
L5 2F 








75 


66 IF 
4i OF 








76 


S5 F 
50 76L 








77 


26 s6 








78 


40 (a * )F 

40 13F 
L3 11S5 


by 44,101 






79 


36 96L 
50 12F 








80 


75 10F 
4o OF 








8l 


L7 OF 
L2 4F 








82 


36 158L 
L5 OF 








83 
84 


66 4f 
S5 F 
32 84l 




Compute unstandardized regression 
weight 




85 


40 (b.jF 
L5 16F 


by 40,102 






86 


74 (x.)F 

1 

40 OF 


by 38,103 






87 


LL OF 
36 88l 




Compute constant term 




88 


26 158L 
L5 OF 
L4 15F 




Change scale 




89 


40 15F 
S5 F 










K 14 



LOCATION 


ORDER 


NOTES PAGE 27 


90 


ko l6F 
50 13F 






91 


75 10F 

kO OF 






92 


L7 OF 
L2 4F 






93 


32 163L 
L5 OF 




Change scale 


94 


66 kF 




Standard error of unstandardized re- 




S5 F 




gression weight 


95 


32 95L 




Waste 


96 


40 (S ) 

F5 48L 
42 48L 


by 41,104 




97 


F5 49L 
42 49L 






98 


F5 51L 
42 50L 






99 


F5 50L 
42 51L 






100 


F5 57L 
42 57L 






101 


F5 77L 
42 77L 






102 


F5 84l 
42 84l 






103 


F5 85L 
42 85L 






io4 


F5 95L 
42 95L 






105 


F5 17F 
40 17F 






106 


LO 6S5 
32 48l 







K 14 



LOCATION 



ORDER 



NOTES 



PAGE 28 



K l*i 



107 

108 
109 
110 
111 

112 

113 
114 

115 
116 

117 
118 

119 
120 
121 
122 

123 
124 



kl OF 
L5 9F 
LO IIP 
50 12S5 
66 9F 
S5 F 
kO 3F 
50 110L 
26 S6 
kO 5F 

L5 5S5 
L4 3F 
kO IF 
LI IF 
32 U.kL 
50 UhL 
26 s6 
50 2F 
7J 10F 
40 4f 
92 139F 
L3 HS5 
32 121L 
L5 4F 
52 ( )F 
50 H9L 

26 SK 
92 131F 
92 515F 
L3 10S5 
32 12kL 
L5 5F 
52 ( )F 
50 123L 
26 SK 
92 131F 



Multiple correlation 
coefficient 



Standard error of estimate 



Waste 



by 7 



Print standard error of 
estimate 



from 118 



by 3 



Print multiple correlation 



LOCATION 


ORDER - 




NOTES 


PAGE 29 


125 


4l 3F 
92 131F 


from 143 






126 


92 515F 
L3 HS5 








127 


32 132L 
L5 ( )F 


by 40, 138 


Print b. . 
-J 




128 


52 ( )F 
50 128L 


by 35, 161 






129 


26 SK 










L5 ( )F 


by 42, 139 ■ 


Print S, 




130 


52 ( )F 
50 130L • 


by 34, 165 


ij 




131 


26 SK 
92 965F 








• 132 


32 132L 
L3 10S5 


by 127 






133 


32 137L 




Print b* . 
ij 






L5 ( )F 


by 43,140 




134 


52 ( )F 
50 134L 


by 33 






135 . 


26 SK 




Print S, * 






L5 ( )F 


by 44, l4l 


ij 




136 


52 ( )F 
50 136L 


by 32, 165 






137 


26 SK 
F5 127L 








138 


42 127L 
F5 129L 








139 


42 12 9L 
F5 133L 








l40 


42 133L 
F5 135L 








i4i 


42 135L • 
F5 3F 








142 


40 3F 
LO 6S5 






_____ _— — — _ — — — — — — — - 



K 14 



LOCATION 



lk} 
lkk 
Ik} 
146 

148 
1^9 
150 

151 
152 
153 
154 

155 
156 

157 
158 

159 



ORDER 



NOTES 



32 125L 
L3 11S5 
32 l^TL 
92 131F 
92 515F 
LI 15F 
52 ( )F 

50 i46l 

26 SK 

92 135F 

F5 25L 

k2 23L 

F5 23L 

k2 25L 

F5 24L 

k2 2kh 

F5 30L 

Fk 16S5 

kO 30L 

L5 14S5 

Lk 14S5 

kO 14S5 

L3 14S5 

32 ( )F 

22 22L 

00 25F 

00 F 

00 1000 0000 OOOOJ 

00 F 

00 0100 0000 OOOOJ 



PAGE 30 



K lk 



by 35,160 



from lkk 
from 25 



Print constant term 



by 



50 7F 
7J 156L 
kO 7F 
lb lk6L 



from 55,82,87 



LOCATION 


ORDER 


NOTES 


PAGE 31 ] 


160 


k2 lk6h 








161 


k6 128L 
L5 13^L 








162 


iA kS5 
k6 l^^L 








163 


26 36L 
L5 130L 


from 73,93 






16k 


iA ^S5 
46 130L 








165 


L5 136L 
iA ^S5 








loo 


k6 136L 
50 8f 








167 


7J 157L 
ko 8f 








168 


26 36L 
00 50SN 

00 850K 











JO 200F 
50 L 








1 


26 S9 

00 421SS 




Store subroutine G on drum 




2 


00 169F 

L5 OF 




[This is a interlude which 
later overwritten.] 


is 


3 


kO kh 
26 999F 

26 850N 








Master Re 


utine 
00 18K 











50 200F 
50 L 




Read subroutine A 




1 


26 S9 
00 SS 




from drum 





K Ik 



LOCATION 



3 
k 

5 
6 

7 
8 

9 

10 
11 
12 

13 
14 

15 
16 

17 
18 



ORDER 



00 kkF 
50 2L 
26 200F 
22 kL 

50 200F 
50 4L 
26 S9 
00 700SS 
00 191F 
22 7L 
50 Sk 
50 7L 
26 S9 
00 900SS 
00 39F 
22 10L 
50 S6 
50 10L 
26 S9 
00 l^OSS 
00 9F 
50 12L 
26 200F 
50 13L 
26 259F 
22 15L 
50 SK 
50 15L 
26 S9 
00 150SS 
00 28F 
22 18L 
50 200F 
50 18L 



NOTES 



PAGE 32 



K Ik 



Enter subroutine A 



Read subroutine B and C and H 
from drum 



Read N-12 from drum 



Read R-l from drum 



Enter subroutine B 



Enter subroutine C 



Read P-l from drum 



Read subroutine D 
from drum 



LOCATION 



ORDER 



NOTES 



PAGE 33 



K Ik 



19 
20 
21 
22 

23 
2k 

25 
26 

27 

28 

29 

30 

31 

32 

33 

3^ 

35 



26 S9 
00 179SS 
00 60F 
50 20L 
26 200F 
22 22L 
50 200F 
50 22L 
26 S9 
00 240SS 
00 29F 
50 2kL 
26 200F 
22 26L 
50 116F 
50 26L 
26 S9 
00 270SS 
00 150F 
50 28L 
26 116F 

22 30L 
50 S6 
50 30L 
26 S9 
00 l^-OSS 
00 9F 
22 33L 
50 SK 
50 33L 
26 S9 
00 150SS 
00 28F 
22 36L 



Enter subroutine D 



Read subroutine E 
from drum 



Enter subroutine E 



Read subroutine F 
from drum 



Enter subroutine F 



Read R-l from drum 



Read P-l from drum 





LOCATION 


ORDER 


NOTES PAGE 3 1 *- K 




36 

37 
38 
39 

List of 



50 200F 
50 36L 
26 S9 
00 teiss 
00 169F 

50 38L 
26 200F 
2^ L 

Constants 

00 58K 

00 (n)F 
00 (n)F 




Read subroutine G 
from drum 

Enter subroutine G 

Stop; Begin program again. 




1 


00 S3 

00 S3 




location of observation row 






00 F 
00 F 




location of means 




3 


00 F 
00 F 




location of sum of cross products 




h 


00 IF 
00 IF 








5 


80 F 
00 F 








6 


80 F 
00 (n)F 








7 


00 F 
00 S7 








8 


00 F 
00 ( )F 




Print parameter for correlations 




9 


00 F 
00 ( )F 




Print parameter for co variances 




10 


00 F 
00 ( )F 




Print parameter for standardized 
regression weights 




11 


00 F 
00 F 




Print parameter for unstandardized 
regression weights 



LOCATION 



ORDER 



NOTES 



PA£E 35 



K Ik 



12 

13 

Ik 

15 
16 

17 



00 F 
00 IF 
00 F 
00 10F 
00 F 
00 F 
00 F 
00 (s)F 
00 F 
00 (n)F 
00 S8 
00 S8 



Dependent variables 



00 ll6K 
I 

Library Routine R-l Square Root 

00 125K 
Library Routine P-l (modified) 





00 


960K 


Library 


Rout 


ine N-12 




00 


86OK 





JO 


S6 




50 L 


1 


26 


S9 




00 


l^OSS 


2 


00 


9F 




L5 


OF 


3 


ko 


12L 




00 


IF 


k 


JO 


SK 




50 


kh 


5 


26 


S9 




00 


150SS 



Print One Number Fractional or 
Integer in a Manner Determined by 
a Program Parameter 



Input a Sequence of Decimal Fractions 



Store R-l on drum 



[This is an interlude which is later 
overwritten . ] 



Store P-l on Drum 



LOCATION 


ORDER 


NOTES PAGE 36 


6 


00 28F 
L5 OF 






7 


^0 13L 
00 IF 






8 


JO sk 








50 8L 




Store N-12 on Drum 


9 


26 S9 
00 900SS 






10 


00 0039F 
L5 OF 






11 


ko i4l 
26 999F 
26 860N 
00 880K 






Library Routine X-7 Sum Check Change FF order in X-7 here to 






read FF 17F 




2k 18N 




— , , _ , — . ., , 



K Ik 



